Azimuthal anisotropy at RHIC: the first and fourth harmonics 
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We report the first observations of the first harmonic (directed flow, vi), and the fourtli harmonic 
(114), in the azimutlial distribution of particles with respect to the reaction plane in Au+Au collisions 
at the Relativistic Heavy Ion Collider (RHIC). Both measurements were done taking advantage of 
the large elliptic flow (1)2) generated at RHIC. From the correlation of 112 with vi it is determined 
that V2 is positive, or in-plane. The integrated V4 is about a factor of 10 smaller than V2. For the 
sixth (ve) and eighth (ug) harmonics upper limits on the magnitudes are reported. 

PACS numbers: 25.75.Ld 



Anisotropic flow, an anisotropy of the particle az- 
imuthal distribution in momentum space with respect to 
the reaction plane, is a sensitive tool in the quest for the 
quark-gluon plasma and the understanding of bulk prop- 
erties of the system created in ultrarelativistic nuclear 
collisions Q. It is commonly studied by measuring the 



Fourier harmonics (vn) of this distribution Elliptic 
flow, V2; is well studied at RHIC and is thought 

to reflect conditions from the early time of the collision. 
Directed flow, Vi, was discovered almost 20 years ago 
and has been extensively studied and reviewed at lower 
beam energies 0- At RHIC energies directed flow in 
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the central rapidity region reflects important features of 
the system evolution from its initial conditions, vi is 
predicted to be small near midrapidity with almost no 
dependence on pseudorapidity. However, it could exhibit 
a characteristic "wiggle" 0, depending on the baryon 
stopping and production mechanisms as well as strong 
space-momentum correlations in the system's evolution. 
A similar rapidity dependence of directed flow could de- 
velop due to a change in the matter compressibility if 
a quark-gluon plasma is formed j^, ^3 ■ I* results in the 
so-called third flow component or "anti flow" ^3 com- 
ponent in the expansion of the matter. This expansion 
direction is opposite to the normal directed flow, vi has 
not previously been reported at RHIC. 

The importance of the higher harmonics in under- 
standing the initial configuration and the system evo- 
lution has been emphasized j^. Recently, Kolb re- 
ported that the magnitude and even the sign of U4 are 
more sensitive than V2 to initial conditions in the hy- 
drodynamic calculations. Those higher harmonics reflect 
the details of the initial conflguration geometry. Besides 
one early measurement at the AGS , reports of higher 
harmonics have not previously been published. 

Experiment — The data come from the reaction 
Au + Au at ^siS" = 200 GeV. The STAR detector [i| 
main time projection chamber (TPC |0|) and two for- 
ward TPCs (FTPC were used in the analysis. For 
the higher harmonics 2 million events in the main TPC 
were analyzed. For the first harmonic analysis there were 
70 thousand events available which included the FTPCs. 

In this analysis the main TPC covered pseudorapidity 
(7?) from -1.2 to 1.2, while two FTPCs covered -4.2 to 
-2.4 and 2.4 to 4.2. The low transverse momentum {pt) 
cutoff was 0.15 GeV/c. In the present work all charged 
particles were analyzed, regardless of their particle type. 
The centrality definition in this paper is the same as used 
previously by STAR [l3|. The errors presented in the 
figures arc statistical. 

Analysis — The difficulties in studying directed fiow 
are that the signal is small and the non-fiow contribution 
to the two-particle azimuthal correlations can be com- 
parable or even larger than the correlations due to flow. 
To suppress the non-flow effects the current analysis uses 
the knowledge about the reaction plane derived from the 
large elliptic flow. One method for eliminating the non- 
flow contribution in a case when the reaction plane is 
known was proposed in 0. It was noted that while the 
correlations of the components of the (first harmonic) 
flow vectors in the reaction plane contain both flow and 
non-flow contributions, the correlations of the compo- 
nents perpendicular to the reaction plane contain only 
non-flow contributions. Then the difference yields the 
flow contribution. Correlating the azimuthal angles of 
two particles {(f>a,4'b), and using the event plane deter- 



mined by elliptic flow {'^2) one gets: 

(cos((/)a - ^'2) C0s((/)6 - ^-2) (1) 
- sin((/)a - ^2) sin(06 - ^z)) 

= (cos(0a + 0h - 2^-2)) ~ Vi^aVl,b{cOs{2{^2 - I'fl.p))), 

where ^'ijp is the azimuthal angle of the reaction plane. 
If only one particle is used to determine the second har- 
monic event plane this expression reduces to 

(cos(0a +(f>b- ~ Vi^aVl.bV2,c, (2) 

which is the basic formula of the three-particle correlation 
method of Borghini, Dinh, and OUitrault The anal- 
ysis of directed flow in this paper is performed using this 
three-particle cumulant method [l8j. The analyses for 
f4, wg, and vg were done relative to the second harmonic 
event plane using the method described in Refs. 0, 01, 
with the event plane resolution calculated from Ref. [3 
equation 11 with fc = 2, 3, or 4. Note that this approach 
in many aspects is very similar to the analysis of directed 
flow described above as it also involves three (for V4, and 
four for vq) particle correlations. For example, for the 
fourth harmonic flow, (approximately, for the exact rela- 
tions actually used in the analysis, see 0) 

(cos(4(/) - 4*2)) ~ vjvi N/2, (3) 

where N is the total number of particles used to deter- 
mine the second harmonic event plane. This expression 
should be compared to Eq. (0). Results obtained with 
this method we designate by V4{EP2}. The anal ysis for 
V4 was also done with three-particle cumulants |2flj by 
measuring (cos(20a + 2(f)b — 40c))- 

vi results — Fig. ^ shows the results in comparison to 
the lower beam energy data of NA49 [13 . The NA49 data 
are also replotted so as to be at the same distance from 
beam rapidity jlil as the STAR results. The RHIC wi(ry) 
results differ greatly from the unshifted SPS data in that 
they are flat near midrapidity and only become signifl- 
cant at the highest rapidities measured. However, when 
plotted in the projectile frame relative to their respective 
beam rapidities, they look similar. It should be noted 
that at the SPS energies of iOA GeV and 158 A GeV HH, 
this y - ybeam scaling does not work, but y/ybeam scal- 
ing does. In the pseudorapidity region \r]\ < 1.2, vi{r]) is 
approximately flat with a slope of (— 0.25± 0.27{stat))% 
per unit of pseudorapidity, which is consistent with pre- 
dictions |i,liE3- 

Note that the sign of vi is undetermined because vi en- 
ters as the square in Eq. (|2Jl. We have plotted vi in the 
positive hemisphere going negative toward beam rapidity 
as it does at the lower beam energy. In the NA49 analy- 
sis j2lj the sign of vi had been determined by defining vi 
for protons near beam rapidity to be positive for periph- 
eral collisions. On the other hand, since the measured 
correlation of Eq. Q is positive, we can conclude that we 
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FIG. 1: (color online). The values of vi (stars) for charged 
particles for 10% to 70% centrality plotted as a function of 
pseudorapidity. Also shown are the results from NA49 (tri- 
angles) for pions from 158^ GeV Pb + Pb midcentral (12.5% 
to 33.5%) collisions plotted as a function of rapidity. The 
open points have been reflected about midrapidity. The NA49 
points have also been shifted (circles) plus or minus by the 
difference in the beam rapidities of the two accelerators. The 
dashed lines indicate midrapidity and RHIC beam rapidity. 
Both results are from analyses involving three-particle cumu- 
lants, 111 {3}. 

have measured the sign of V2 to be positive. While the 
absolute values of V2 at RHIC arc well determined 
this is the first direct indication that the elliptic flow at 
RHIC is in-plane. 




1 2 3 4 5 6 

p, (GeV/c) 

FIG. 2: (color online). The minimum bias values of V2, V4, 
and V6 with respect to the second harmonic event plane as a 
function of pt for |?7|<1.2. The V2 values have been divided by 
a factor of two to fit on scale. Also shown are the three particle 
cumulant values (triangles) for V4 (1)4(3}). The dashed curves 
are 1.2 • V2 and 1.2 • uf. 



Vi results — The results as a function of pt are shown 
in Fig.|5]for minimum bias collisions (0 — 80% centrality). 
Shown for V4 are both the analysis relative to the second 
harmonic event plane, V4{EP2}, and the three-particle 
cumulant, W4{3}. Both methods determine the sign of V4 
to be positive. As a function of pt, Vi rises more slowly 
from the origin than but does flatten out at high pt 
like V2- The ve{pt) values are consistent with zero. The 
hydrodynamic calculations of Kolb ^3 for pions from b = 
7 fm collisions agree very well with our measured V4 for 
charged particles for centrality 20 to 30%. However, he 
calculates to be —1.2% at 2 GeV/c, while we observe in 
Fig.[21for minimum bias data that it is essentially zero. It 
also appears to be zero in our data for all the individual 
centralities. Ollitrault has proposed for the higher 

harmonics that Vn might be proportional to ^2^^ if the 
(j) distribution is a smooth, slowly varying function of 
cos(2(/)). In order to test the applicability of this scaling 
we have also plotted V2 and t;| in the figure as dashed 
lines. The proportionality constant has been taken to be 
1.2 in order to fit the V4 data. 

Kolb ^3 points out that for 112 > 10%, which occurs 
at high Pt, and no other harmonics, the azimuthal distri- 
bution is not elliptic, but becomes "peanut" shaped. He 
calculates the amount of V4 (which looks like a four-leaf 
clover) needed to eliminate this waist. Our values of V4 
as a function of pt arc about a factor of two larger than 
needed to just eliminate this waist. 

The results for V4 as a function of pseudorapidity are 
approximately flat in the acceptance of the main TPC 
(|ry|<1.2) with an average value of (0.44 ± 0.02)%. How- 
ever, in the FTPCs (2.7<|r;|<4.0) the average value is 
(0.06 ± 0.07)%, consistent with zero, with a two sigma 
upper limit of 0.2%. Consistent with the first observa- 
tion by PHOBOS [3, at 77 = 3 for minimum bias colli- 
sions we observe V2 = (3.06 ± 0.10)%, which is a factor 
1.8 smaller than at midrapidity. Thus V4 seems to fall off 
faster at high rapidity than V2 ■ This faster fall off at high 
pseudorapidity is also consistent with V4 scaling like w|. 

Fig. 131 shows the centrality dependence for pt- 
integrated W2, V4, and vq with respect to the second har- 
monic event plane and also V4 from three-particle cumu- 
lants (i;4{3}). The five-particle cumulant, W4{5}, (not 
shown in the figure) is consistent with both methods but 
the error bars are about two times larger. The vq val- 
ues are close to zero for all centralities. These results 
are averaged over pt, thus refiecting mainly the low pt 
region where the yield is large, and also averaged over rj 
for the midrapidity region accessible to the STAR TPC 
(|77|<1.2). To again test the applicability of Wj^^ scal- 
ing we have also plotted and the figure as dot- 
ted histograms. The proportionality constant has been 
taken to be 1.4 to approximately fit the V4 data. The 
larger constant here compared to that used in Fig. |21 is 
understood as coming from the use of the square of the 
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FIG. 3: (color online). The pt- and 77- integrated values of 
V2, V4, and vq as a function of centrality. The 112 values have 
been divided by a factor of four to fit on scale. Also shown are 
the three particle cumulant values for V4 (1)4(3}). The dotted 
histograms are lA ■ V2 and 1.4 • u|. 

average instead of the average of the square, and because 
the integrated values yield-weight low pt more, where the 
best factor is slightly larger. 

The Vn{EP2} values averaged over pt and 77 (|77|<1.2), 
and also centrality (minimum bias, — 80%), are (in 
percent) V2 = 5.18 ± 0.005, V4, = 0.44 ± 0.009, vq = 
0.043 ± 0.037, and = -0.06 ± 0.14. Since is es- 
sentially zero, we place a two sigma upper limit on wg of 
0.1%. Also, iig is zero, but the error is larger because the 
sensitivity decreases as the harmonic order increases. 

Systematic uncertainties — In both approaches, f4{3} 
and Vi{EP2}, the non-flow effects are suppressed com- 
pared to the case where the fourth harmonic event plane 
is used. The remaining non-flow correlations, along 
with event-by-event flow fluctuations, are thought to be 
the major contributors to the systematic uncertainties. 
Background from secondary particles is expected to be 
less than 15%, and remaining acceptance effects are mea- 
sured to be very small. All errors and limits quoted so far 
are statistical, and should be increased by the systematic 
uncertainties below. 

From non-flow effects we estimate the relative system- 
atic uncertainty in W4{3} to be about 20%. The largest 
contribution comes from situations in Eq. (jSJ where one 
particle is correlated with one of the other particles due 
to non-flow, and with the third particle via flow. Our 
estimate is based on the assumption that the entire dif- 
ference in the published values |3| of V2{EP2} and W2{4} 
is due to non-flow effects. Comparison of W4{3} to f4{5} 
leads to a similar estimate for this systematic error. 

Prom non-flow effects we estimate the relative system- 
atic uncertainty in wi{3} also to be about 20%. Our es- 



timate is based on the assumption that our two-particle 
correlation value of vi using only the first harmonic event 
plane in the FTPCs, vi{EPi}, of about 3% is entirely due 
to non-flow effects. 

The other effect important for the comparison of our 
results to theoretical calculations is event-by-event flow 
fluctuations. As was discussed Q, flow measurements are 
done by two or many particle correlations, resulting in, 

not {vn), but (v^y^'^ . If flow fluctuates event-by-event, it 
could lead to a difference between these two quantities. 
Fluctuations in the initial geometry of the collision at 
fixed impact parameter can account for the difference be- 
tween V2{Ep2} and W2{4} 0, and also between v^lEP^} 
and W4{3} [23j. Although the flow fluctuation contribu- 
tion to W4{3} is greatly reduced, it still could lead to an 
effect of about a factor of 1.2 to 1.5. 

Conclusions — We have presented the first measure- 
ment of vi at RHIC energies. ^1(77) is found to be approx- 
imately flat in the midrapidity region, which is consistent 
with microscopic transport models, as well as hydrody- 
namical models where the flatness is associated with the 
development of the expansion in the direction opposite 
to the normal directed flow. Within errors we do not 
observe a wiggle in vi (ry) at midrapidity. The pseudora- 
pidity dependence of vi in the projectile fragmentation 
region is very similar to that observed at full SPS energy. 
We observe a positive correlation between the first and 
second harmonics, indicating that elliptic flow is in-plane. 
This is the first direct measurement at RHIC of the ori- 
entation of elliptic flow relative to the reaction plane. 

We have measured V4 as a function of pt, 77, and cen- 
trality. We observe that appears to scale approxi- 
mately as as a function of pt, 77, and centrality. wg, 
although essentially zero, is not inconsistent with scaling 
as f 2 . This is the first measurement of higher harmonics 
at RHIC and it is expected that these higher harmon- 
ics will be a sensitive test of the initial configuration of 
the system, since they provide a Fourier analysis of the 
shape in momentum space which can be related back to 
the initial shape in configuration space. In fact, it has 
been emphasized that V4 has a stronger potential than 
V2 to constrain model calculations and carries valuable 
information on the dynamical evolution of the system. 
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